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Abstract We study the nature of fluctuations in variety of price indices involving companies listed on the New York 
Stock Exchange. The fluctuations at multiple scales are extracted through the use of wavelets belonging to Daubechies 
basis. The fact that these basis sets satisfy vanishing moments conditions makes them ideal to extract local polynomial 
trends, through the low pass or 'average coefficients'. Subtracting the trends from the original time series yields the 
fluctuations, at different scales, depending on the level of low -pass coefficients used for finding the 'average behav- 
ior' . The fluctuations are then studied using wavelet based multifractal detrended fluctuation analysis to analyze their 
self-similar and non-statistical properties. Due to the multifractality of such time series, they deviate from Gaussian 
behavior in different frequency regimes. Their departure from random matrix theory predictions in such regimes is also 
analyzed. These deviations and non-statistical properties of the fluctuations can be instrumental in throwing significant 
light on the dynamics of financial markets. 



1 Introduction 

Financial time-series which were in the past of interest to only economists, have led to considerable inter-disciplinary 
research due to the applicability of various physical laws and techniques in their analysis. This has led to the discovery 
of various new aspects like fractality (1), multifractality (2), correlated behavior (3, 4) and complex network structure 
(5). The last few years has also seen a lot of activity in terms of explaining the correlations in financial markets through 
the Random Matrix Theory (RMT) framework (6-8). 

Fractals as first predicted by Benoit Mandelbrot in the 1960s (9) have been widely applied to understand various 
processes in Physics (10), Chemistry (11) and Biology (12). Mandelbrot in 1963 proposed the study of fluctuations 
in the market prices (1) which opened new vistas for the analysis of stock markets through statistical physics. In 
the recent years, Peng et al. proposed the "Detrended Fluctuation Analysis (DFA)" in 1994 (2) to study the DNA 
nucleotide structure as a random walk problem which was extended to study the price fluctuations in economic time- 
series under a mono-fractal hypothesis (13-15). However, the inadequacy of the mono-fractal hypothesis to model 
the behavior of financial time series was soon pointed out and consequently, a multi-fractal model called the "Multi- 
Fractal Detrended Fluctuation Analysis (MFDFA)" was proposed (16). This method used a variable window approach 
to calculate the local variances in the profile of the data series from the polynomial fit. Manimaran et al. in 2009(17) 
building on the work by Kantelhardt et al. proposed the "Wavelet Based Multi-Fractal De-trended Fluctuation Analysis 
(WB MFDFA)", where, using the Multi-Resolution Analysis (MRA) capable "fractal like" kernels, the time-frequency 
resolution and extraction of fluctuations for multi-fractal analysis was shown to have a greater efficiency than its 
predecessor. 

The study of the correlation matrix of the financial return series have been shown to agree well with the predictions 
of RMT and the nearest-neighbor-spacing of the rank-ordered unfolded eigenvalues of the correlation matrix follow 
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that of the Gaussian Orthogonal Ensemble (GOE) (7, 18). This behavior dubbed as an universal behavior of the 
financial return series have again, recently been studied with respect to temporal evolution of financial correlations to 
study the differences between the assumption of strongly correlated financial time-series and uncorrelated financial 
time-series (19). In this context, it becomes important to analyze the nature of correlations in the time-frequency 
domain in order to ascertain the effects of non-stationarity and transience on such studies. The inefficient handling 
of such signals by either Fourier Transform or Short Term Fourier Transform for the purpose of a time-frequency 
localized study have already been established, leading to the development of Wavelet Transform (20). 

In this work, we will use a wavelet based fluctuation extraction technique to study the correlations of the fluctuations 
at various frequency windows (called scales in the wavelet parlance). We will also briefly comment on the multi-fractal 
nature of the fluctuations and the distributions of the associated parameters: Hurst exponent and singularity strength. 
The organization of the article is as follows: in Sec. (2), we will briefly review the theoretical methods of Wavelet Based 
Fluctuation Extraction (WBFE), Wavelet Based Multi-Fractal De-trended Fluctuation Analysis (WBMFDFA) and the 
Random Matrix Theory (RMT) based method to study the time-frequency localized correlations of the fluctuations. 
Further in Sec. (3), we discuss the results obtained by the application of WBFE, WBMFDFA and correlation analysis 
on the price index of 196 scrips trading on the New York Stock Exchange (NYSE) between September 1984 to June 
2010. Finally, we summarize and conclude with the scope for future work in Sec. (4). 



2 Review of theoretical methods 

2.1 Wavelet Based Fluctuation Analysis 

In the following, the analysis of a time series given by X(t) is carried out by calculating the "log-normalized return 
series" R(t): 

r(f)=lDgX(»+l)-logX(t) (1) 
R(t) = r{t) - {r{t)) (2) 

where, (•} and a r are the time average and standard deviation of the log-return series r(t). a r is also called as the 
"volatility of returns". The profile Y(t) is calculated by taking the cumulative sum of the log-normalized return series: 

r(t) = £*(*) (3) 
*=i 

which is then subjected to the Wavelet Based Fluctuation Extraction (WBFE)(17). The WBFE can be performed 
following the steps: 

1. Calculate the one-dimensional discrete wavelet transform (lDWT)(20-23) of the profile: 

/ 

Y(?)= £ + £ L d abWab{t) (4) 

b=-oo a>0b=-°° 

where q, are the "low-pass" coefficients that capture the trend or the average behavior of the signal and d u \, are the 
"high-pass" coefficients capturing the local fluctuations in the signal at various window sizes a. The functions <P 
and f are called the "scaling filter" and the "high pass" filters respectively. The father and mother wavelets <j>t,(t) 
and y/ a b(t) are orthogonal to each other and are subjected to the admissibility conditions 



J <$>{t)dt < oo (5) 

J y{t)dt = (6) 

I fy*(t)y(t)dt = 0, orthogonality (7) 

J \<j)(t)\ 2 dt= 1 = J \\j/{t)\ 2 dt (8) 
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The scaled and translated versions of y/(t) are called the "daughter wavelets" 

Yabit) = 2" /2 \i/{2 a t-b),a eR,beZ+ (9) 

which differ from the mother wavelet \jf(t) at the a th scale by 2 a in height and 2"l 2 in width, a and b are called 
the scaling and translation parameters respectively and I = [logA^/log2j is the maximum number of scales for the 
profile Y(t) of length N. 

The wavelet kernel for the 1DWT should be chosen such that it captures the maximum information from the signal. 
For example, the Daubechies' family of wavelets satisfy vanishing moment conditions which make them blind 
to various polynomial trends. The wavelet Db-N (with the index number N being even integers between 2 and 
20), has N/2 vanishing moments limiting the representation of a polynomial trend of N/2 in the signal. The Db-4 
wavelet has two vanishing moments making it blind to constant and linear trends. In this work the Db-4 wavelet is 
employed. 

2. Calculate the approximate trend T a (t) at the scale of interest a and subtract it from the profile Y(t) to get the 
fluctuations Z„(f): 

Z a {t)=Y(t)-T a {t). (10) 

The Z a (t ) obtained by this method represent the fluctuations at different frequency bands (the scale is inversely related 
to the frequency). Consequently, these fluctuations can be probed to analyze the behavior of the signal in various 
frequency bands like Fourier power law and moments of the fluctuation distribution. It has been shown earlier that the 
well-known /~ 3 behavior of market fluctuations appear only in the low frequency or long wavelength regime(24). Due 
to the varying window sizes (corresponding to different scales) and the convolution error generic to wavelet transforms, 
the extracted fluctuations happen to have erroneous values at the edges. These errors are corrected by performing the 
WBFE on the reversed profile and then taking the average of the two (forward and reversed) fluctuation series at each 
scale. 



2.2 Wavelet Based Multi-Fractal De-trended Fluctuation Analysis 

These fluctuations can also be subjected to a multi-fractal analysis which is a modified form of the original MFDFA 
proposed by Kantelhardt et al. in 2002(16). 

The fluctuations obtained in Eq. (10) are further subdivided intoA( s = \N /s\ segments of size s such that s = 2 a ~ l W, 
where W is the support width of the wavelet and a represents the scale. Thus the fluctuations obtained at various scales 
can be analyzed at window sizes corresponding to the scale and the wavelet used. 

Since the fluctuations are guaranteed to have zero mean, we can directly find the variance of each segment and thus 
calculate the fluctuation function 



F g (s) 



q^0 



(11) 



F q (s) is the "gth order" fluctuation function, where q e [—m,m], m 6 Z + . The negative (positive) q values capture the 
fractality of the broader (finer) fluctuations. It can be easily seen that at q — 0, the Fq(s) blows up. Hence, to circumvent 
this issue, at q = 0, 



F (s) = exp 



— 2>g{d 2 ( M )} 



?/2 



k=\ 



1/'/ 



q = 0. 



(12) 



The "generalized Hurst exponent" h(q) can be obtained from F q (s) since, as a function of the segment size s, F q (s) 



follows a power law of the form 



Mq) 



(13) 



It must be noted that at q = 2, this method reduces to the standard mono-fractal fluctuation analysis and h(2) is the 
Hurst exponent. The multi-fractal scaling exponent i(q) can be calculated as 



x(q)=qh(q)-l 



(14) 



The singularity spectrum f(j5) is related to the multi-fractal scaling exponent t(q) by a Legendre transform 
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1-Z{q), and / ( ±z(q)) = q^x{q) - %{q) = f{fi) = qfi x{q) (15) 
dq \dq J dq 



2.3 Correlation analysis of fluctuations 

The fluctuations Z a obtained at scale a through Eq. (10) can be analyzed for correlations. The fluctuations for the entire 
data set of N scrips can be written as a JV x I matrix JT ', where T is the time length of the fluctuation series of each 
the N scrips. The correlation matrix is thus given by 



¥{x,y) 



1 



Sty Sty" 1 



(16) 



where, • • - T is the transposition operation. In the case where SC consisted of N mutually independent normally dis- 
tributed fluctuations of length T, c £ could be considered to be a Wishart matrix (18, 25, 26). Under the constraint 
N — > ooj T — > °° and Q = T/N> 1 the density of eigenvalues of the correlation matrix takes the form (27) 



P(A) 



Q V / ( A max-A)(A -A min ) 



A 



with a\ : being the variance of the matrix 3£ and A^ x given by 



(17) 



A, 



max _2 

- a or 



Q 



■±2* 



(18) 



It must be noted that under the constraint — > °°, the eigenvalues of the the matrix ^ lie strictly in the range 
[A m i n ,A m in]. However, for finite sized matrices, there exists a finite probability of finding eigenvalues outside this 
range. Indeed, it has been shown that for financial time series, the nearest-neighbor eigenvalue spacing A = A, +J t — A,, 
obtained through unfolding the eigenvalues of ^ follow the distribution for a Gaussian Orthogonal Ensemble (GOE)(7) 



(19) 



Here, we will investigate the scale dependence of the p(X) for the correlation of the fluctuations Z„(f) and Z' a (t) at 
scale a, where Z a (t) and Z' a (t) are different scrips in the corpus. This will give us an idea about the effectiveness of the 
RMT predictions in different frequency regimes. 



3 Results and discussion 



3.1 Data 



We have analyzed 196 scrips trading on the New York Stock Exchange (hereafter referred to as NYSE) in the period 
from September the 7th, 1984 to June the 10th, 2010. The scrips analyzed in this work have been chosen so as to 
encompass high-cap, mid-cap as well as low-cap sections of the American stock market. This combined with the time- 
frame for analysis includes some of the major crashes of the NYSE, for the example the Black Monday on October 
ls9, 1987, the July 2, 1997 crash triggered by the Asian financial crisis, the burst of the dot-com bubble leading to a 
three years sluggish activities from March 10, 2000 and the two years long bear run of the market from 2007-2009. 
The data has been analyzed using the WBMFDFA method and also through a correlation analysis of the fluctuations 
in the random matrix theory framework. 
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Fig. 1 Plots depicting the density of unfolded eigenvalues p(A) of the correlation matrix of the fluctuations at (a) different scales and at 
(b) scale 10. It can be seen the p(A) fits with the analytical Eq. (17). 



3.2 Time-frequency localized correlation analysis 



Analysis and investigation of correlations of the fluctuations over different frequency windows a{aoc \ Jf where / is 
the frequency) can provide us with insight into the spectral behavior of the market correlations thereby improving our 
understanding of the collective behavior of the market in different time spans. For example, if the correlations between 
different scrips representing different sectors of the market in short time windows (high frequency, low scale) is low, 
then this could be exploited to guard the simultaneous crashes of different sectors in the event of a crisis. However, 
if the scrips are correlated in the low frequency limit, then the long term correlations of the market could lead to the 
"healing" of the market in a systematic and efficient way after a crash. This kind of information could be very useful 
for the policy makers in order to both monitor the economy as well provide for safe-guards against possible crashes. 



200,, 




(a) (b) 

Fig. 2 Plots depicting the density of nearest neighbor unfolded eigenvalue spacing p(A) of the correlation matrix of the fluctuations at (a) 
different scales and at (b) scale 5, ps(A) fit with p(A) = alexp (— foA 2 ) where a = 35.64 and b = 0.7707 (with 95% confidence bounds) 
as an illustration. We can see that though the GOE fits well in the high and mid frequency ranges, at low frequency or long periods, they do 
not fit with the RMT assumptions. 
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In order to understand the nature of spectral correlations in the market, we analyze the spectrum of the correlation 
matrices at different scales. In Fig. (1), we have shown the density of the unfolded eigenvalues of the correlation 
matrices of the fluctuations obtained at different scales. We observe from Fig. (1(a)) that at lower and middle scales 
(upto a = 7), the unfolded eigenvalue distribution p(A) is largely uniform and changes to fit the Eq. (17) at higher 
scales as shown in Fig. (1(b)). Since in this analysis, we have looked at the behavior of the fluctuations over the 
whole time period T = 5799, the value of Q = 29.58 is very large. We could expect that the correlations of such 
fluctuations will show significant deviations from the RMT prediction since it fits well in short time windows(28). We 
must remember that the fluctuations obtained by the WBFE method at characteristic of the frequency range (scale) 
under study. The spectral correlations under investigation here show that at high frequencies (low scales), the GOE fit 
of Eq. (17) is not followed by p(A). However, at higher scales, that is at lower frequencies, the p(A) fits reasonably 
well with the Eq. (17). 

In Fig. 2(a), we have plotted the density of the eigenvalue spacing p (A) against different scales. It is clearly visible 
that at the lower and middle scales, the p(A) follows the RMT prediction well, which we have exemplified in Fig. 
2(b) at scale 5. This is in contrast with the results in Fig. (1(a)) suggesting that the eigenvalue spacing at lower scales 
is behaves differently from the eigenvalues themselves. This is an interesting observation which we believe should be 
explored in greater details. 



4 Conclusion 

To summarize, we have analyzed the nature of fluctuations from different sectors of the New York Stock Exchange 
(NYSE) through Wavelet based multi-fractal analysis and RMT based techniques. The number of companies being 
small, it is expected that there will be significant deviation from the RMT prediction. Interestingly, although the 
density of unfolded eigenvalues exhibited this behavior, the nearest-neighbor eigenvalue spacing distribution matched 
reasonably well with the RMT predictions at lower and middle scales. This points out that the spacing distribution can 
shed considerable light on the nature of the high frequency fluctuation correlation, for a smaller corpus of data. 

In conclusion, wavelet analysis when combined with RMT approach can reveal considerable information about the 
correlations at different scales and is quite useful for modeling the behavior of high and low frequency components of 
physical processes. 
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